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ABSTRACT 


A mean streamline analysis of the effect of stator blade orienta- 
tion on the performance of an axial flow compressor was performed by 
means of a computer program. Measurements were made on a 3-stage axial 
flow compressor at the Naval Postgraduate School at six stator stagger 
angles between 23.8° and 44.3° for a fixed orientation of the rotor 
blades. Maximum efficiency and pressure ratio were measured at a 
stator stagger angle of 31.8°. Results at other blade settings showed 
that by varying stator stagger angle with flow rate optimum efficien- 
cies and pressure ratios can be achieved over a wide range of operating 
conditions. 

The results of the analysis were compared with the measured results. 
Suggestions are made for improving the manner of adapting cascade test 
data to performance predictions. 

By applying a non-dimensional deflection coefficient it could be 
shown that minimum work input corresponded to maximum efficiency. 

The test compressor has a tip diameter of 36 in. and a hub/tip 
ratio of 0.6. The blading tested is of the free-vortex type with a 


design degree of reaction of 0.5. Tip speed was about 185 ft/sec. 
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SECTION 1 


INTRODUCTION 


Axial flow compressors have a wide range of application. They 
are usually matched to a turbine in a set. Both open-cycle sets, 
such as aircraft turbo-jet engines, and closed-cycle sets for power 
generation are in service. In every application compressor efficiency 
is the critical factor in set performance. Accurate predictions of 
off-design compressor performance must be made during the design stage 
to be able to make corrections prior to manufacturing. 

Theoretical analysis of the flow in an axial compressor involves 
three-dimensional partial differential equations of a complex nature (1). 
In these equations it is difficult to account for real gas effects such 
as boundary layer growth on machine walls and blade surfaces. The 
resulting wakes behind blade rows create non-uniform conditions; hence, 
it becomes necessary to base prediction methods on experimental data. 

These data would be obtained best on an actual compressor. How- 
ever, most compressors have relatively short blades and small flow 
annuli. Even if only small pressure probes were inserted between the 
rows of blades, the flow in these machines would be substantially 
altered. Test rigs must therefore be used that have large dimensions. 

A first approximation of the flow in axial turbomachines can be 
obtained in rectilinear cascade test rigs. The intersections of a 
co-axial cylinder with the blades of a row produce a series of iden- 
tical and identically oriented profiles which are unwrapped into a 
plane to establish the corresponding rectilinear cascade. A finite 


number of straight blades with profiles similar to, but larger than 
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those obtained,is then arranged in a rectangular duct through which 

air is blown at the appropriate inflow angle. Flow surveys taken ahead 
of and behind the blade row give experimental data of air turning angle 
and total pressure loss for various blade geometries. Such a recti- 
linear cascade is in use at the Naval Postgraduate School. The results 
of extensive testing of compressor airfoil shapes at NASA have been 
summarized by Lieblein (2). 

The usual problem in compressor design consists in selecting flow 
areas, blade shapes and blade layout to satisfy design specifications. 
Additional studies are frequently necessary if the design performance 
is not reached. An example is the study of Vavra which was made when 
the so-called Clark CSN-1l compressor of the ML-1 nuclear gas turbine 
of Aerojet-General Nucleonics failed to perform as required. Vavra 
analyzed the design changes proposed by the manufacturers and made 
further recommendations (3). However the ML-1l project was cancelled 
before the improved compressor could be built so that it was not 
possible to verify the suggested changes by experiments. 

The so-called inverse problem consists in predicting the perform- 
ance of an existing machine with a known blading by means of rectilinear 
cascade data and other experience factors. Such an approach was carried 
out by Gibbons and Bartels (4) for the 12-stage Allis-Chalmers axial 
flow compressor which supplies air to the turbine test facilities at 
the Naval Postgraduate School. Because of the unorthodox performance 
of the first stages, and also due to small flow annuli and short blade 
heights, it was not possible to predict the performance of the machine 
accurately. However, the analysis gave important indications for 


design changes to improve the performance of the compressor. 
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This thesis makes another attempt to solve the inverse problem 
for a 3-stage axial flow compressor built by the California Institute 
of Technology with the support of the Office of Naval Research. The 
compressor is sufficiently large so that the insertion of pressure 
probes causes relatively small flow perturbations. The machine was 
designed to permit variations in blade shapes, blade angles, tip 
clearances, and staging. The installation at the California Institute 
of Technology is described by Bowen, et al. (5) This report also 
describes a theory of perfect fluid flow in axial flow turbomachines. 
The results of this theory are compared with the measured data of the 
first stage of a so-called "free-vortex" type blading. Part 2 of the 


report by Bowen, et al. (6), gave a detailed investigation of multi- 


stage flow for both "free-vortex" blading and a more complex "solid 
body rotation" blading. Measurements of blade skin friction losses 
were greater than expected from cascade tests. The growth of boundary 
layer along the inner and the outer annulus walls of the flow was less 
than expected. Alsworth and Iura (7) conducted extensive and detailed 
measurements of flow patterns in a single stage of 'free-vortex" 
blading. They carried out accurate measurements of the blade skin 
friction losses and the radial distribution of work input. The final 
report from the California Institute of Technology was a hot-wire 
anemometer study of compressor stall by Iura and Rannie (8). They 
observed that stalled flow regions rotate in the direction of blade 
rotation without changing shape but with a speed that is not propor- 
tional to rotor speed. 

After the Turbo-Propulsion Laboratory of the Naval Postgraduate 


School was built, the compressor was relocated there. Since then 


Z3 


it has been used for laboratory courses to supplement instruction 
in basic theories of turbomachines. 

The present study was conducted for this compressor because of the 
possibility of changing its blade angles. For fixed stator and rotor 
blade angles a highly peaked curve of efficiency versus flow rate is 
usually obtained with axial flow compressors. At a slightly changed 
stator blade angle the peak will be displaced somewhat. A family of 
efficiency curves for various stator stagger angles is expected to 
yield a flatter efficiency curve over a wider range of flow rates for 
fixed rotor blade angles. The wide variation in flow rates required 
of the engines for the supersonic transport has led to design pro- 
posals by General Electric for compressors where the stator blade 
angles can be changed during operation of the jet engine. Hence, the 
thesis topic is of current interest. 

For the tests it was necessary to obtain extremely accurate 
pressure measurements. An extremely sensitive bourdon tube pressure 
indicator was added to the instrumentation. A permanent Pitot-static 
tube was installed, and accurate flow rate calibrations were carried 
out. A parallel theoretical effort produced predictions of com- 
pressor performance by an existing computer program which was modified 
and adapted for use on the I. B. M. Model 360 computer system installed 
at the Naval Postgraduate School. 

The author wishes to express a debt of gratitude to the faculty 
of the Department of Aeronautics of the Naval Postgraduate School for 
his engineering education. Particular thanks are due to Professor 
M. H. Vavra for his enthusiasm and guidance in showing what can be 


done with that education. 
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SECTION 2 


O. N. R. 3-STAGE AXIAL FLOW COMPRESSOR 


The compressor installation is shown in Fig. 1. The compressor 
inlet is in the background of the picture. The exit section is 
equipped with a throttle valve. A torque meter is installed on the 
shaft between the drive motor and the compressor. The overall dimen- 
sions of the equipment are given in Fig. 2. 

A basic criterion of the design of the compressor was flexibility 
of operation. Each blade row may be removed. The angle setting of 
each blade is adjustable in 0.5° increments. In order to minimize 
flow disturbances by pressure probes the outside diameter of the blad- 
ing is 36 in. For the same reason, the compressor has a hub/tip ratio 
of 0.6, giving an inside diameter of 21.6 in. and a blade height of 
7.2 in. To permit simulation of conditions in a multistage unit, the 
machine has three stages. One row of inlet guide vanes simulates the 
effects of previous stages. Two rows of exit guide vanes remove the 
whirl component after the third stage. To insure that general align- 
ment and tip clearances are maintained the compressor casing is very 
rigid. 

The outer casing consists of two cast iron half-cylindrical 
shells bolted together along the horizontal plane through the axis. 
The rows of inlet guide vanes, stator blades, and exit guide vanes 
are held by bolts extending through the casing (Fig. 3). Inside the 
casing a hollow steel shaft carries three cast iron drums for each 
of the rows of rotor blades. The shaft is supported by two roller 


bearings whose outer races are pressed into the supporting strut 
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assemblies (Fig. 4). These assemblies rest in close fitting channels 
in the casing (Fig. 5). The six assembly struts have symmetrical air- 
foil sections. Figure 6 shows the nine blade rows after assembly. 

The casing was designed to give maximum accessibility for mea- 
suring instruments. Six rectangular instrument ports are located in 
the upper half at 30° from either side of the vertical (Fig. 3). They 
are placed in the second and third rotor planes, the first, second and 
third stator planes, and behind the third stator. The ports hold a 
special instrument carriage permitting detailed flow surveys in any 
axial plane. The casing has numerous radial survey holes whose loca- 
tions are specified in Fig. 7. 

The blades are made of ALCOA 356 aluminum alloy without heat treat- 
ment. There are thirty rotor blades and thirty-two stator blades per 
row. The details of construction of the blades are illustrated in 
Figs. 8 and 9. The blading is designed to have a degree of reaction 
of 0.5 at the mean radius. The rotor blades are twisted by 49° from 
hub to tip, and the stator blades by oe At the mean radius, the 


important blade parameters are: 


Parameter Rotor Stator 
re) fe) 
Camber angle 20. 32 29.98 
Design stagger angle 43.80° 28.80° 
Thickness-to-chord ratio 0.1 0.1 
Chord length 2.60 in. 2.60 in. 


The values of the blade parameters at other radii, and the method 
used for calculating the thickness distribution, are specified in 


Ref. 5. 
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The point of maximum thickness was set at 0.35 of the chord 
length from the leading edge. The thickness was modified over the 
rear 15 per cent of the section to provide a trailing edge thickness 
0.02 in. The thickness distribution is applied about a parabolic 
mean camber line. Figure 10 is a picture of a "free-vortex" stator 
blade. Tip clearances of 0.020 in. for stator, and 0.037 in. for 
rotor, are maintained by the use of shims at the point of attachment. 

The attaching device permits variation of rotor and stator blade 
Stagger angle. Figure 11 shows how these changes are made. The blade 
shaft extends through the rotor drum at the hub. An adjusting plate 
is secured to the blade shaft by a tapered pin. The plate has a series 
of eleven taps on a circle about the shaft axis which are spaced at 
intervals of 4.5°. Between the adjusting plate and the rotor drum is 
a fixed sector which is attached to the drum by a pin which guarantees 
its alignment. The fixed sector has eleven holes which are spaced at 
intervals of 4.0° of arc about the blade shaft axis. Blade alignment 
is maintained by a set screw through both holes in the adjusting plate 
and the fixed sector. If the center holes of the plate and the sector 
are lined up, the blade is at the design stagger angle of 43.8° for 
the "free-vortex" rotor. An angle change of 4.0° can be accomplished 
by keeping the set screw in the center hole of the adjusting plate and 
moving the blade so the screw is inserted into the next hole of the 
fixed sector. Small angle changes of 0.5° are accomplished by moving 
the blade slightly so that the two holes immediately adjacent to the 
center holes are lined up since blade angle is changed only by the 
difference in arc between the two holes. A similar arrangement is used 
for the stator blade settings (Fig. 10). Figure 12 illustrates stator 


blade settings which are one degree smaller than the design stagger 


ay 


angle. The two blades on the right-hand side of the figure belong to 
the first stator row. The set screw is two holes away from the center 
or reference holes. The blades on the left-hand side of Fig. 12 belong 
to the row of inlet guide vanes. They are set at an angle which is by 
0.5° smaller than the design stagger angle. 

The assessories to the basic compressor will be described by fol- 
lowing a path from the inlet to the exit and to the power source. The 
instrumentation of the compressor will be described in its appropriate 
place in the same sequence. Details of instrument calibrations are 
contained in Appendix A. 

The inlet duct is seen in the background of Fig. 1. It consists 
of a screen, an entrance bellmouth, and a length of straight pipe. 
Figure 13 shows the large mesh screen which prevents the ingestion of 
foreign matter from the compressor bay apron. Mounted on the screen 
is a mercury thermometer used to determine ambient temperature. [It 
has provisions for psychrometric analysis of the incoming air. 
Immediately behind the screen is the bellmouth. It changes the dia- 
meter rather abruptly from almost two diameters at the flare, to 
36 in. in the axial distance of about 15 in. A cylindrical inlet duct 
two diameters long connects the bellmouth to the compressor proper. 

The inner surface of the duct is enameled to give smooth flow surfaces. 

The flow rate through the compressor was determined by surveys in 
the inlet duct. Details of the calibration are given in Appendix B. 
The survey plane is about midway between the bellmouth and the com- 
pressor. The surveys were made with a Prandtl-type Pitot tube. [In 
Fig. 14 the probe is installed for taking a vertical traverse. It is 


possible to disassemble the probe while the compressor is running (Fig. 15). 
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It may then be remounted for a horizontal traverse. Traverses at 30° 
and 60° from the vertical direction are possible. These traverses can 
be made by unbolting the entire inlet duct from the compressor and 
rotating it on its cradle. The design of the probe prevents measure- 
ments closer than 0.25 in. to the inlet duct wall. 

A honeycomb straightener is installed behind the inlet Pitot tube 
to equalize the flow entering the compressor (Fig. 14). Attached to 
the forward strut assembly is an ogival wooden fairing (Fig. 5) to 
provide a smooth transition of the flow from the inlet duct to the 
compressor annulus. 

Pressure measurements may be taken at any one of the radial 
survey ports described. Figure 16 shows the probe used for these 
measurements. It is a United Sensor three-hole probe, model YC-120. 

A central hole measures stagnation or total pressure. A static port is 
located on each of two faces of a wedge at an angle of 45° from the 
total pressure hole. The static pressures are balanced on a water 
manometer to insure alignment of the probe in flow direction. The 
probe is mounted on a probe holder which has vernier scales permitting 
radius adjustments to 0.01 in. and angle adjustments to 0.1°. 

The accurate determination of flow rate is essential to compressor 
analysis. For this purpose a modified Prandtl type Pitot-static probe 
was installed. Figure 17 shows the installation details. The probe 
is inserted in a radial survey hole located ahead of the inlet guide 
vane row. It protrudes forward into a channel between two adjacent 
struts of the supporting assembly which is shown in Fig. 5. The probe 
is approximately at the center of the channel between struts. The 


probe is positioned in the radial survey hole by a brass plug, and 
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secured to the casing by a simple locking device. Figure 18 shows the 
probe and its locking device in the background. In the middle of the 
figure are seen the three-hole probe and its holder located between the 
inlet guide vanes and the first rotor. In the foreground, a plug, 
originally located at the same axial position as the probe, has been 
removed to show a hole through which the probe can be inserted. 

Pressures obtained by the probes were measured by a Texas Instru- 
ments Fused Quartz Precision Pressure Gauge which is shown in Fig. 19. 
The velocity head in the inlet is about 0.6 in. of water. A low- 
pressure bourdon tube was obtained for the gauge which measures from 
O to 166 in. of water with the instrument reading from 0 to 200,000 
counts. The calibration constant of the low pressure bourdon tube was 
determined to be 240.423 counts/psf of the pressure gauge (See Appendix 
A.). The pressure gauge operates both in manual and servo modes, the 
latter providing automatic nulling of the unit. 

Since flexibility of pressure selection was desirable, two 
Giannini Sp-101lA pressure scanners were connected to the pressure and 
reference ports of the bourdon tube. One of tle scanners is seen in 
the foreground of Fig. 19. The arrangement of pressures to the 12- 
channel switches is shown in Fig. 20. Later in the study a simpler 
system of five valves and two manifolds was constructed (Fig. 21) to 
eliminate leakage flows that seem to have occurred in the scanners. 

At this time the original 3-hole probe was placed permanently behind 
the third stator (Fig. 22). Another 3-hole probe was placed ahead of 
the first stage rotor. This arrangement permitted direct measurement 


of the pressure increase across the three stages ,A Pao: 
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Figure 22 also shows the instrument traverse carriage mentioned 
before. It was used for a survey of total pressure in peripheral 
direction behind the third stator at two different radii (See Appendix 
C.). It is possible to vary the radial locations of the probe with an 
accuracy of 0.01 in. The probe may be rotated through 360° about its 
axis. The carriage can move the probe 15° in peripheral direction, 
which covers a whole blade spacing since a stator blade channel covers 
11.25°, and a rotor blade channel 12.0° in peripheral direction. The 
locations of the radial survey holes with reference to the indicated 
meridional angle of the carriage are shown in Fig. 23. The represented 
Stator profiles are at the design stagger angle at the mean radius. 
The actual stagger angle measured with respect to a plane through the 
axis is the complement of the indicated angle. 

A short cylindrical duct connects the compressor to the exit 
elbow,which is equipped with sheet metal turning vanes to minimize 
losses. A transition piece connects the elbow to the throttle valve 
which is seen in Fig. 1. The valve consists of two rectangular metal 
doors which move in slides. The position of the doors is controlled by 
a right- and left-hand threaded lead screw which is rotated by a 
Graham variable-speed transmission. A revolution counter indicates 
the approximate valve opening. 

The compressor is driven by a 50 HP Fairbanks Morse induction 
motor. It requires a three-phase,sixty-cycle, 440 volt power supply, 
and can operate at two fixed speeds; namely, at about 900 rpm or 
1200 rpm. 

Attached to the motor shaft is a Baldwin-Lima-Hamilton SR-4 


torque meter type A. It is shown in Fig. 24 during calibration with 


I 


known weights attached to a lever. In the background may be seen a 
Brown Instruments strain gauge readout which is of the standard 
Wheatstone bridge type. A conversion constant of 4.1565 ft-1lb/1b of 
torque meter readout was determined (See Appendix A.). Under the 
protective cover next to the torque meter in Fig. 24 is installed a 
six lobe flux cutter to measure the motor: speed by means of an 


electronic counter. 
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SECTION 3 


FLOW RATE CALIBRATION 


The permanent Pitot-static tube was installed to provide rapid 
and accurate measurements of volume flow rate. At a given throttle 
valve setting this probe was calibrated against the result of two 
traverses of the cylindrical inlet duct. Through-flow velocities in 
the flow annulus ahead of the inlet guide vanes were varied from 
110 ft/sec to 70 ft/sec, so that compressibility effects could be 
ignored. The objective of the calibration was to determine a cali- 
bration constant which gave the volume flow rate when multiplied by 
the velocity measured at the permanent probe. 

An initial survey of the inlet duct was made at the highest pos- 
sible flow rate of 512 fie egal Figure 25 shows a plot of the mea- 
sured velocities of both horizontal and vertical traverses. The 
smooth velocity profile in the duct is evident. Figure 26 is a graph 
of the measured velocities in the inlet duct at the lower flow rate of 
310 PP lieecs which shows a higher degree of scatter. For these tests 
the pressure indicator was operated in the servo mode. Even in the 
meter mode the indicated pressure oscillated considerably. Pressure 
fluctuations seemed to have been initiated by the vibration of the 
brass rod on which the traversing Pitot-static tube was mounted in the 
inlet duct. The pressure indicator corrects a disparity between the 
angular position of a mirror of the bourdon tube and the indicator 
dial at the rate of the full scale reading of 200,000 counts in 120 
seconds. With the oscillating pressure applied to the sensitive 


instrument, oscillations of 50 counts were observed for a pressure 


a2 


corresponding to a reading of 600 counts. To reduce these oscillations, 
variable length capillary tube damping devices were installed both in 

the total and static pressure lines of the traverse Pitot-static tube. 
Details of the arrangement and the results of tests are shown in 
Appendix B-l. 

A special survey was made to determine the thickness of the bound- 
ary layer on the walls of the cylindrical inlet duct. However, con- 
struction of the traverse probe prevented taking readings closer than 
0.25 in. from the wall. Figure 27 is a plot of the readings at high 
and low flow rates. In each case it is evident that large losses occur 
near the walls. With the velocity distributions of Fig. 27 the flow 
rate is about 2.5 per cent smaller than the value obtained without 
considering the changes near the walls. 

Data were taken in the boundary layer at intervals of 0.25 in. 
from 17.75 in. to 16.00 in. In the main stream,data were taken at 
4.0 in. intervals from 16.00 in. to the centerline. With horizontal 
and vertical traverses, four points were obtained at each radius and 
two at the axis for a total of 46 locations. At each station the 


following measurements were taken: 


q. - velocity head in the inlet duct (counts) 
qo - velocity head of the permanent Pitot-static probe (counts) 
Poi - static pressure with reference to atmosphere in inlet 


duct (counts) 
r - static pressure with reference to atmosphere at permanent 


Pitot-static probe (counts) 


~ - temperature at the inlet screen (°F) 
Pay - barometric pressure (in. of Hg) 
Tear - temperature of column of mercury (°F) 
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The barometric pressure P r (in. Hg) was converted into pounds 


Ba 


per square foot to obtain the atmospheric pressure PA by 


fas (71. 467) (ps fa) (1) 


Bar 
where the constant was obtained by correcting the height of the column 
of mercury for temperature variations in specific gravity of the liquid. 
The inlet temperature was converted to absolute temperature,and was 
considered to be the total temperature Te because of the low veloc- 
ities, or 

ie, sealers +> 207 (PR) (2) 
The static pressure P. at the inlet was converted from the meter 


reading to an absolute pressure by 


P.; 
P. a ae (psfa) (3) 
Fn 240,423 


where the constant was obtained by the calibration procedure described 

in Appendix A. The pressure Py measured at the permanent probe was 

obtained in a similar manner. The local air density a at the inlet 
u 


was computed from 


~ On00ASTO GE DNC Fe (6 gexio” (ib pec” /ét”) (4) 
(4.696144) Teo 9 Lg 


and similarly Po at the permanent probe. The velocity head in the 


inlet q,, was converted from the meter reading, 


ae 
~ pyeepgeeewencn £ (5) 
bis 240.423 a.) 


The dynamic head at the permanent probe q,, was obtained similarly. 


The velocity V5 in the inlet is then 


Vio =7y _ (£t/sec) (6) 
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A similar procedure gave the velocity Nee at the permanent probe. The 
arithmetic average of the 46 values of Wi obtained during the carrying 


out of the traverses gives the average value Ne of 
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| 
wa ( iro (ft/sec) 
a 


The measured values of Vv; were adjusted to account for variations in 


casey = 


flow rate due to slight compressor speed variations 


Ni ge 
NN, =a p= - (ft/sec) 
© 


The four values of V obtained at each radius were averaged by 


iad 
4. 
a | 
Maia 4 as Viad} (ft/sec) 
re 
The volume flow rate was first calculated by neglecting the losses in 


the boundary layer to give V+ - In general the volume flow rate 


is obtained by 


F = 27 V Ir dr Gemeiece) 


oO 


Lf the above equation were integrated directly, the velocity at the 
axis would not be included in the volume flow rate calculation 
because the zero radius is a member of the product under the integral. 


Therefore the integral was replaced by a summation. 


— « 3 

Nae 2 ), Viay aA (£t~/sec) 

The quantity OA = 2 1 r h/144 represented that part of the annulus 
area in the center of which V ia was measured, where h is the height 


of the area. This method assumed that the velocity measured at the 


radius of 16.0 in. exists also at the wall. <A volume flow rate 


36 


deficiency was then calculated to account for the difference in 
velocities at 16.0 in. radius and those measured in the boundary 
layer, or 

Vp ® Viav nt 7 Viav et. (£t/sec) 
By a trapezoidal integration the volume flow deficiency Vp becomes 


18,0 


van a 2a Vp ia dr ieee 
16.0 


Then the actual volume flow rate is 


ae wir a ve (ee eee) 


Se 
A calibration constant CCC for each run was obtained from 


CCC :———— (£t/sec per ft/sec) 
Oav 

A first series of calibration runs was invalidated when a leakage 
hole was discovered in the traversing Pitot-static tube. After repairs, 
the calibration constants CCC were obtained with a maximum relative 
error of 0.1 per cent. The final calibration constant was 4.4050 
a per ft/sec velocity obtained from the readings of the perma- 
nently installed Pitot-static probe. 

A data reduction computer program, ONRFLO, for these calibration 
procedures,was generated for the I. B. M. Model 360 computer. A 


listing of the program is included as Appendix B-2. Sample output 


data from the final run is given in Appendix B-3. 
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SECTION 4 


MEASUREMENT OF COMPRESSOR PERFORMANCE 


The objective of the study was to determine the effect of changes 
of the stator blade orientation on compressor performance. [It was 
decided to restrict performance investigation to the three stages only; 
namely, from a station ahead of the first rotor to a station after the 
last stator, or between locations SP-2 and SP-8 of Fig. 7. The effects 
of the inlet, the inlet guide vanes, the exit guide vanes, and the 
diffuser were not considered. 

For tests at a particular stator stagger angle the barometric 
pressure Paar was determined first. For various settings of the 
throttle valve the following measurements were taken: 

T - ambient temperature (°F) 


re) 

P so - static pressure with reference to atmosphere at the 
permanent Pitot-static probe (counts) 

q - velocity head at the permanent Pitot-static probe (counts) 


P - total pressure with reference to atmospheric between the 


inlet guide vanes and the first rotor (counts) 


ras P,_,~ total pressure difference across the three stages (counts) 
TY - reading of torque meter (lb) 
N - compressor speed measured by electronic counter (rpm) 


With the relations listed on p. 35 these measuring data establish 


the following quantities: 


Py = Poe (71.467) (ps fa) (1) 

a = 1. + 459.7 (°R) (2) 
Pso 

P= Rea (ps fa) (3) 


38 


p, = 01002376 Te SSL. Peis t2x \ hee) | 


(14,696)( 144) ae ae, 


‘bo 

bor ~ 340,423 a, i 
2. 

Vo = e. : (ft/sec) (6) 


The absolute total pressure ahead of the first stage, or after the 


inlet guide vane, Been is 
=. ae (ps£) (7) 
. s 
tiga ‘ 240, 423 
The pressure difference across the three stages A P 3 ca is obtained 
from 
aP 3.5 
Ae : a (psfa) 
3-3A 240.423 


From the torque meter reading T, the actual torque TY. is 


Vyrq Neem ees) ) (£t-1bf) 


where the constant was obtained by the calibration procedure described 


in Appendix A. The angular velocity W is with the measured speed, N 


27 N 


MS moms (rad/sec) 


Both volume and weight flow rates were calculated. The volume 


flow rate Vee is then 


3 
Ve oma Vin (O44 cea) (ft/sec) (8) 
where the constant was obtained by the calibration procedure described 


in the preceding section and Appendix B. The weight flow rate w is 
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32.174 Fa 518.7 - & 
Ve cana = Ve 7, (7.886) (1b0/sec) re 


For comparison with the prediction program a referred flow rate we is 


defined by 
° W Liga 
Wi = Jee ae (1bm V°R/sec psia) (10) 


FP. /144 
Two basic indices of performance were computed; namely, the 


3-stage pressure ratio re , where 


3-s 


Pers s + 4 en 


K-3-§ = —————— (11) 
Feiga 
and the 3-stage total-to-total efficiency “i obtained from 
on 
iss 2 ——— (12) 
Pwa 
where Poth and Pa are the theoretical and actual horsepower, respec- 
tively. The theoretical power is Poth even by 
ior = Ww ah ;, (778, 3) (ft-1bf/sec) 
where Ah. . is the isentropic specific enthalpy difference for a 
frictionless compression from P.. to P_, after the third stator. 
tiga td 


Assuming a perfect gas there is 


a 


ahi. = .90Sip al, 3-5 (Btu/1bm, °R) 
where AT 4. is the isentropic temperature difference from the total 


ahead of the first stage, to T / after the third 


temperature 
P BLennn td 


iga 
stator (See Fig. 28.). Since Ceies very nearly equals the ambient 


temperature Tho? there is 
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T y 
Ah, ~ ae Ie a - ] (Btu/1bm, °R) 
tiga 





With the isentropic relation 


/ fal 
le a | ae ¥ 
Neiga Keiga 


the final expression for the determination of the theoretical power is 


| x! 
Far, = W alia )( et ]778.3 (£t-1b£/sec) (13) 


The actual power is obtained from 
— = Tre Ud (ft-lbf/sec) 


The total-to-total efficiency at the three stages is then 


Yok rv 
eases Sete | (Res) ‘— } 778, 3 (14) 


Wg 
For comparison with other compressor data certain non-dimensional 
parameters are introduced. The average flow function Pay is the 


ratio of through flow velocity Ny and peripheral speed at the mean 


diameter U__, or with 
av 


Nar 


Ri Nepali} amr. aS 


and, for the mean diameter of 14.4 in. 
14.4 
Va” Gare 
aN 12, (ft/sec) 


there is 


12)(144) Fi 
P ay = os. ——— (0.1842) (15) 
OD (144)(651.4) OO 
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An average dimensionless stage pressure coefficient itn, is defined by 


Ve 
AT ae = YW CE ie (psf) 


or 





Ww = er 


AV 


2 (16) 
3 es av 


The so-called dimensionless deflection coefficient be is defined by 


APs sth _ OP ssa _ ty po US ast) 
3 > 2a-5 "es 


Tove Uw (pst) 


} 


Hence, 


v _ Wor 


av ~~ el 
13-5 


The performance calculations were made for the conditions along 


(17) 


the mean radius of the stage,14.4 in., by assuming that these conditions 
are representative of the performance of the 3-stage compressor. 
Detailed investigations which support this simplification are given in 
Appendix C. The data reduction was accomplished by the computer 
program ONRETA, which is listed also in Appendix C. 

Data were taken on 11 occasions for the following settings of 


stator stagger angle: 


Run Stator Stagger Angle 
1,2 7 a 
4, 5 he al 
6, 7 31.8 
Bens 35.8° 
10 39.8° 
ll 44.3° 
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SECTION 5 


PERFORMANCE PREDICTION PROGRAM 


The program used to predict the performance of multi-stage com- 
pressors is called AXCO3. It is listed in Appendix D and represents 
an adaptation of a program by Vavra (3), modified by Gibbons and 
Bartels (4). 

The method assumes that the flow is axisymmetric. Computations 
are made for the streamlines at the mean radius of a stage. A stage 
by stage analysis of the machine is performed. The output of one 
stage is used as input for the next stage. The computation is started 
at a flow rate less than that anticipated for surge, and the flow rate 
is then increased by small increments to cover the whole operating 
range of the machine. 

Program AXCO3 uses two principal subroutines; namely, one called 
ROTOR and the other called STATOR. They establish the velocity 
triangles of a rotor or a stator row of blades, respectively, from the 
known geometry of the bladings, flow rate, and discharge conditions of 
the preceding row of blades. Both subroutines make use of three 
additional subroutines. 

Subroutine THEORY calcualates the flow conditions for minimum 
profile losses. Cascade data by Lieblein are used for this purpose (2). 
These data are presented as curves for zero-camber incidence angle 
Gus) , slope factor n, zero camber deviation angle (€,) Io? and 
slope factor m, all given as functions of the inlet air angles 6 , 


or 1? with the blading solidity O° as a parameter. The solidity 


O is defined as the ratio of blade chord c and blade spacing s. 
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The curves described above have been expressed by polynomials in powers 
of solidity and inlet air angle which are used in subroutine THEORY. 
The accuracy of the polynomials has been verified, and the errors are 
listed in Appendix D. The incidence angle for minimum loss is deter- 
mined by an iteration procedure which changes the inlet air angle until 
agreement is reached within a tolerance of 0.02°. The deflection which 
occurs at the minimum profile loss can then be calculated with the 
established polynomials also. 

The subordinate subroutine CASCAD calculates the performance of 
the blading at the actual incidence or inlet flow angle of the blade 
row by uSing a graph of Lieblein that establishes the change of the 
deviation angle with incidence angle (d 5 /d .) as a funettonvot 
inlet flow angle and solidity. 

In the earlier computer programs the curves of this graph were 
expressed by a polynomial that covered the range of solidities from 
0.0 to 1.8. Since the curves have nearly exponential character, which 
cannot be expressed with ease by polynomial functions, it was found 
that an error of the order of 10 per cent was possible. Therefore, 

a new polynomial was established for solidities between 0.6 and 1.4 
which limits the error to less than 1.9 per cent. 

Subroutine CASCAD further utilizes experimental limits on blading 
performance by Howell (9) which are illustrated in Fig. 29. Howell 
found that for most cascades the ratio of the profile drag coefficient 
Cy at an arbitrary incidence angle and Comin at the minimum loss 
incidence a is a unique function of the quantity Guar e., ; 
where € io is the optimum flow deflection for the incidence angle Lee 


Since the applicability of the data presented in Ref. 2 is limited to 
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design point calculations, where the incidence angles are close to the 
nominal incidence ae the curves cannot be used for off-design analy- 
ses without restrictions. Without such limits the flow deflections € 
could be increased simply by increasing the incidence angle ters since 
the data of Ref. 2 do not establish criteria for flow separations 

which are associated with a radjcal increase in profile losses. Exper- 
imental results of Howell show that the flow deflection € at incidence 
L should not exceed the optimum deflection €,, by more than 25 per cent 
to avoid excessive losses. He states that this limit occurs if the 
profile loss coefficient Ch is about twice the minimum loss coefficient 
Comin’ The solid curve in Fig. 29, labeled Ee /eé 190» is supposed to be 
the relation between E/E 10 and (Os LigPME 1 as obtained from subroutine 
CASCAD with the data from Ref. 2. It will be assumed that this curve 
holds for negative values of (.L= Lay Gay, only. For positive values 
of this parameter the dashed curve will be applied to satisfy the 
limitations of Howell. This curve is somewhat arbitrary but represents 
the best estimate for the actual operating condition of compressor 
cascades. Figure 29 further shows a curve relating the ratio ColiC ie 
to (Lali € rp » which is adapted from Ref. 9. Both this curve and 
the dashed curve for the deflection ratio have been expressed analy- 
tically in subroutine CASCAD. If the ratio C,/C) ain exceeds a value of 
2.0 for positive quantities Gis Y/Er6 ,» the program prints out 
“surge in blade row xx" and proceeds with calculations at an increased 
flow rate. The positive value of (U- 4/15 at which this occurs is 
the input quantity SMAX. If Ch/C) in exceeds a value of 2.0 at 


incidence angles lL smaller than Cio, the program prints out "minimum 


pressure in blade row xx" and halts computation. The negative value of 
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Gye) /e (5 at which this occurs is the input quantity SMI. In 
both instances the symbol xx refers to the number of the blade row, 
and indications are given also whether the row is a stator or a rotor. 
For any incidence angle within the useful operating range,the values 
of the deflection ratio and of the ratio = are computed by sub- 
routine CASCAD and control returned to the blade row subroutines 
ROTOR or STATOR. 

The efficiency of a blade row is determined by subroutine ETACAL. 
Losses due to tip clearances, secondary flows, and wall friction are 
taken into account with the relations proposed by Vavra (1). 

If conditions of surge exist, the flow rate is increased by about 
1.0 per cent and the entire computation is repeated. If the incidence 
angles are in the so-called useful operating region of Fig. 29, the 
output of subroutine ROTOR or STATOR becomes the input for the sub- 
routine that calculates the next row of blades. Interstage data are 
printed out as the program progresses from one stage to the next. If 
all blade rows of a machine have been processed, the overall compressor 
efficiency and pressure ratio are computed for the particular flow rate. 
This flow rate is then increased by a specified amount, and the calcu- 
lating process is repeated. If the so-called minimum pressure is 
reached, the computation stops and a summary of the overall performance 
parameters is printed out. 

The details of the calculation of the flow through the inlet duct, 
the inlet guide vanes, the exit guide vanes, the diffuser and the dis- 
charge passages are not described in this thesis since it is concerned 
only with the performance prediction of the three stages of the com- 


pressor from a station ahead of the first rotor to the discharge at 
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the third stator. The methods applied for the flow analysis in these 
passages may be found in Vavra (3) and Gibbons and Bartels (4). 
According to Howell there occur additional effects that influence 
the off-design calculations. The growth of boundary layers on the 
walls of the annular flow channel will change the velocity profiles 
and can reduce significantly the effective flow area. The program 
takes account of this effect by a blockage factor that can vary from 
stage to stage. This boundary layer growth is responsible also for 
increasing peaks in the radial distribution of the axial velocity 
components in succesive stages. These higher velocities outside of 
the wall boundary layers produce smaller incidence angles for the 
main portion of the flow, thereby decreasing the actual flow deflection 
and reducing the work absorbed by the fluid in the stage. Hence 
in a compressor that consists of stages with identical bladings, 
the last stages will produce smaller pressure ratios than the stages 
at the compressor inlet. Experience shows the efficiency in successive 
stages is not reduced by the peaking of the velocity profiles, and 
their effects are usually taken into account by a so-called work-done 
factor that is about unity for the first stage and gradually decreases 
for the successive stages. Such work-done factors can be introduced 
in program AXCO3, but for the present compressor they will be taken as 
unity because of the small number of stages and the large blade 
heights. For the same reasons the blockage factor is assumed to be 


equal to unity also. 


The geometry of the blades could be obtained from the information 
of Ref. 5. Values at the mean radius were used for the analysis. 


Bowen (5) determined that the average flow angle at the mean radius 
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after the inlet guide vanes ALE was 20.0° for the design inlet guide 
vane stagger angle of ibe independent of flow rate. From Fig. 12 
it can be recognized that the inlet guide vanes were set at a stagger 
angle of 10.9° for the tests. Therefore a value of 19.5° was used in 
the analysis program for the average flow angle leaving the inlet 
guide vanes. Reference 5 does not give values of Comin for the 
blading. It is possible to introduce in the program various values of 
Comin both for rotor and stator blading. Hence it was possible to 
estimate an average apparent value of Comin for blading by comparing 
measured results to predictions of the analysis program for several 
assumed values of Comin’ 

A prediction of the 3-stage efficiency for different stator 
stagger angles is shown in Fig. 30 for an assumed value of Comin of 
0.020. The light line connecting the individual peak efficiencies 
represents the calculated operating envelope which could be obtained 
if continuous control of stator stagger angle for maximum efficiency 
were possible. 

For comparison with measured performance a series of calculations 

O O O 


were made for stator stagger angles of pao L168: 0gp cd Oe ay 


39.8°, and 44.3. For each stagger angle, values of C of 0.000, 


Dmin 
0.006, and 0.008 were assumed. To avoid confusion with the computa- 
tions of actual performance based on measured quantities, which have 
been called runs, each calculation of the analytical prediction 


program has been assigned a two-digit case number. The first digit 


identifies the stator stagger angle of the case as follows: 
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First Digit Stator Stagger Angle 


1 23.8° 
2 2756" 
3 31.8° 
, 35.8° 
5 39.8° 
6 44. 3° 


The second digit identifies the assumed value of Comin of the case 
as follows: 


Assumed Value of C 


Second Digit Dmin 
1 0.000 
2 0.006 
3 0.008 


For example, case 42 is the predicted performance at a stator stagger 
angle of 35.8° with a value of Comin of 0.006. The summary output of 


each case is presented in Appendix D. 
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SECTION 6 


DISCUSSION OF RESULTS 


Experimental data for performance measurements were taken 
during 55 hours of running time. Another 81 hours of operation were 
used for miscellaneous calibration and testing, including 44 hours 
for flow rate calibrations. 

The rotor stagger angle was set at 43.8° for all runs. The 
design stator stagger angle is 28.8. For the performance tests, the 
stator stagger angle was set at 208°; Paap oh 31a. 35.8, 39.8°, 
and 44.3°. Measured performance is presented in comparison with the 
predicted results of AXCO3,and also by establishing dimensionless 
performance parameters. Figures 31 through 36 are graphs of measured 
performance in comparison with the prediction by AXCO3 for each stator 
blade angle setting. Figures 37 and 38 are summary plots of maximum 
efficiency and pressure ratio for all stator blade angle settings. 
Figures 39 through 44 are graphs of the non-dimensional parameters at 
each stator stagger angle. Figure 45 is a summary plot of maximum 
efficiencies and associated deflection coefficients for all stator 
blade angle settings. Figure 46 is a summary plot of maximum pressure 
coefficients for all stator blade angle settings. 

On each of Fig. 31 through Fig. 36 the 3-stage efficiency reaches 
a peak value for each stator stagger angle. The referred flow rate at 
which this peak occurs decreases as the stator stagger angle increases. 

Also plotted on each of Fig. 31 through Fig. 34 are the results 


predicted by AXCO3 for the three estimated values of C of 0.000, 


Dmin 
0.006 and 0.008. The predicted 3-stage efficiency reaches a peak for 
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each Crees. The referred flow rate associated with predicted peak 
efficiency decreases with increase in stator stagger angle. The 
referred flow rate at which the predicted peak efficiency occurs 
decreases with increase in Comin at a given stator stagger angle. 

AXCO3 predicts the condition of surge for each case according 
to the criteria of Howell shown in Fig. 29. The surge condition is 
assumed to exist whenever the parameter Wiig exceeds the 
maximum positive value SMAX which has been inserted in the program. 
The incidence angle for the rotor blade rows becomes larger as the 
through-flow velocity decreases until a surge condition is indicated. 
AXCO3 predicts surge at the same referred flow rate for each setting 
of stator stagger angle. This is due to the fixed angle of the 
inlet guide vanes. The assumption has been made that the air flow 
angle leaving a blade row does not change with variations in flow 
rate in the incompressible flow regime. This result was verified by 
Bowen (5). The inlet guide vanes were not rotated as the stator 
stagger angle was varied. At all stator stagger angle settings the 
first stage rotor had the same incidence angle at particular flow 
rates. In the stage-by-stage analysis, "surge in rotor 1" was 
indicated at the same flow rate in each run, regardless of the stator 
stagger angle setting of subsequent stators. The region of referred 
flow rates smaller than that at which surge is indicated by the 
program is marked with a dashed line on Figs. 31 through 34, and is 
less than 41.5 lay Onveee psia for every case. 

When the value of the parameter (L-Cie-.. is determined to be 
less than SMI, the negative value of the parameter which has been 


inserted in AXCO3, the program stops the computation and prints 
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"minimum pressure in blade row xx." The referred flow rate at which 
computation is stopped decreases as stator stagger angle is increased. 
At stator stagger angles of 31.8° and lower (Figs. 31 through 33) the 
referred flow rate at which computation is stopped is larger than the 
maximum flow rate of the graph. At a stator stagger angle of 35.8° 
(Fig. 34), the computation stops at a referred flow rate of 45.5 

Liem YOR/see psia. At stator stagger angles of 39.8° and 44.3° (Figs. 
35 and 36), the program did not establish useful data. At the high 
stator angles the referred flow rate at which computation was stopped 
was lower than the referred flow rate for surge in the first stage 
rotor. The latter flow rate remained unchanged at a fixed inlet guide 
vane angle. 

This difficulty illustrates some of the limitations of program 
AXCO3 when applied to an actual machine. The program assumes that 
the compressor is designed properly. Changes from stage to stage 
must be gradual. For example the annulus area may change gradually 
from stage to stage with no adverse effects. The checks for minimum 
pressure and surge limits are applied in each blade row. Surge 
in any blade row is interpreted as surge in the machine, and similar 
limitations are imposed for minimum pressure. The program cannot 
cope with discontinuities similar to those that occur by leaving 
the angle of the inlet guide vanes unchanged. It is possible to 
imagine a situation where the angle setting of the inlet guide vanes 
might have to remain unchanged; for instance, if the inlet guide 
vanes must support the front bearing of the rotor shaft. In such a 
design the first rotor would have higher blade efficiencies than 


succeeding blade rows at large flow rates,while at lower flow rates 


a2 


the blade efficiencies of the first rotor would be smaller than for 
succeeding blade rows. The 3-stage efficiency of an actual machine 

is an overall parameter which includes the effect of the different 
blading efficiencies for different blade rows. In an actual machine 
the flow will adjust itself to these conditions,whereas program AXCO3 
cannot cope with these peculiar circumstances. For these reasons the 
program did not produce performance data for the higher stator stagger 
angles of 39.8° and 44.3°. 

The detection of the actual surge point during the tests is 
difficult with the available instrumentation and must be based on 
acoustic phenomena. Incipient surge,probably due to rotating stall, 
was associated with an unmistakable oscillating change in sound level 
produced by the compressor. This effect appeared suddenly even when 
flow rate was decreased slowly. To eliminate the sound due to this 
surge condition it was necessary to increase the flow rate by about 
5 to 10 per cent. A flow rate slightly above the surge point could 
then be reached by gradual throttling. For each run one test point 
was taken in the surge region and one test point as close as possible 
to the surge point. Efficiency at test points in the surge region was 
below 0.90 for all stator stagger angles. The efficiency at the data 
point close to stall rose from 0.89 to 0.94 as the stator stagger 
angle was increased from 23.8° to 31.8°, and efficiency decreased to 
0.89 with further increase in stator stagger angle to 44.3°. As the 
surge condition is approached the rate of decrease of efficiency is 
set by the relative location on the efficiency versus referred flow 
rate graph of the two data points in the surge region and close to 


surge. At a stator stagger angle of 31.8° which produced the maximum 
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3-stage efficiency and at S608F , efficiency decreased most rapidly as 
the surge condition was approached. For other stator stagger angles 
the efficiency decreased less rapidly. It is not possible to make a 
meaningful comparison between measured surge point data and the surge 
condition calculated by AXCO3 because surge was predicted to occur at 
the same referred flow rate for each case. 

It is possible to compare measured values of maximum 3-stage 
efficiency with the maxima calculated by program AXCO3. For all 
four stagger angles at which computations were made, the maximum 
3-stage efficiencies were predicted to occur at flow rates about 
5 per cent greater than the flow rates at which the measured maxima 
occurred. According to Howell the point of maximum efficiency 
occurs for a value of the parameter qi- Ces )/E;q about equal to 
+0.19 (Fig. 29). At greater positive values of the parameter the 
losses are believed to increase rapidly. Since the measured maximum 
efficiency occurred at a flow rate which was about 5 per cent less 
than was predicted, an analysis of off-design performance of a com- 
pressor stage was made. The method of Vavra (Ref. 1) assumes that 
the air flow angles leaving a blade row do not change with flow rate 
in incompressible flow. The velocity triangle changes caused by a 
reduction from design flow rate are shown in Fig. 47. The subscript 
d refers to the design condition, and the primed quantities refer to 
off-design. The analysis was conducted at the maximum efficiency 
point for a stator stagger angle of 27.8. The following quantities 
were measured: 


0.517 


Design flow coefficient Po 


Inlet flow angle A, (fixed) 20.0 
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The off-design flow rate was fixed at 5 per cent less than the design 
flow rate,to correspond to the flow rate at which the measured maximum 


efficiencies occurred. The off-design flow coefficient is 
/ 


ae 
i 


(0.95) (0.517) 


The peripheral components of velocity are: 














V 
o, = q tan x, 
a 
a 
oa es =a ees 
Wout z 
Nut cee = 1- 0.95 d, tan a, 


The relative inlet flow angles for both conditions are: 


sida arctan] — fan a, | 
aeban | walling ~ tan | 


yt 


6,, = arctan 


fa 
ae = arctan Wor /U 


0.95 J 0.95 


d 
With the same rotor stagger angle, the difference in incidence 


angle Ll between design flow rate and off-design flow rate is the 
difference in OF 

ee: Uy ue 6° - C,, Z 59,.2°- aes 5° — ae Oey Ps 
The difference between the parameter CABIN EH at the design 
flow rate and the parameter at a flow rate 5 per cent less than 


design is 


A (t-bo) - falls = +015 
Elo 11.4 


where €,, was calculated by AXCO3. The measured maximum efficiencies 
occurred at a value of the parameter (t-¢)))/€,, which was greater 
by 0.15 than the value of the parameter for which the program AXC0O3 
predicted maximum efficiencies. 

Lieblein's data were obtained from compilations and reductions of 


cascade test data primarily for blade profiles which were NACA 


oe, 


65 (A - series airfoils built up on equivalent circular arc 


10? 
camber lines. The free-vortex blades installed in the compressor 
have a thickness distribution which was derived from the theory of 
Ref. 1 applied to a parabolic mean camber line. If the actual 
incidence for minimum profile loss for the installed blades were 
smaller by 1.7° than calculated by subroutine THEORY, all values of 
the parameter (ae ens calculated in subroutine CASCAD would be 
larger by 0.15 than they actually are. The flow rate at which the 
maximum efficiency is calculated would be 5 per cent greater than 
the flow rate at which maximum efficiency was actually measured. 

Howell presented results of cascade tests of British C-1l and 
C-2 airfoil shapes for several values of solidity. The Ch /C) in 
curve of Fig. 29 was obtained by interpolation for a solidity of 
0.8. If the upper limit of the parameter (UL )/E for minimum 
loss were extended by 0.15 from 0.19 to 0.34, the flow rate at which 
maximum efficiency was calculated analytically would coincide with 
the flow rate for measured maximum efficiency. 

The measured efficiency curves are more peaked than the curve 
computed by AXCO3. If the Ch /Ch an curve of Fig. 29 had a greater 


slope in regions where CH /C approaches a value of 2.0, the pre- 


Dmin 
dicted efficiency curves would be less flat and would more nearly 
approximate the curves of measured performance. 

It was one object of this thesis to attempt to determine the 
value of Cer for the installed blade profile. Vavra (1), p. 377, 
indicates that these values usually do not exceed 0.008. For each 


stator stagger angle, analytical calculations were made for estimated 


Cc. ain of 0.000, 0.006 and 0.008. The results for the smaller stator 
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stagger angles are plotted onFig. 31 through Fig. 34. The value of 


Comin was estimated for each stator stagger angle: 
Stator Stagger Angle BS eee Corus 

23.8 0.008 

21748 0.004 

31.8 0.003 

35.8 0.007 


The arithmetic average of the estimated values is 0.0055, but there is 
is obviously considerable scatter in the experimental data. 

A summary plot. of maximum peak efficiencies is presented in Fig. 
37. If continuous control of stator stagger angle were possible, as in 
the proposed SST engine, Fig. 37 would be an operating curve for max- 
imizing efficiency at any flow rate. The peak maximum 3-stage efficien- 
cy was 0.956, measured at a stator stagger angle of 31.8. Also shown is 
a flatter curve of predicted peak 3-stage efficiencies using Comin of 
0.006 for the four lower stator stagger angles. The maximum predicted 
value is 0.942 at a stator stagger angle of 30°. The predicted values 
occur at flow rates 5 per cent higher than corresponding measured val- 
ues. The measured values decline more steeply at stagger angles other 
than the angles for maximum efficiency than do the analytically calcu- 
lated values. The curve relating stator stagger angle to performance 
shows that,in the region near the maximum efficiency, relatively small 
changes in angle are required for large variations in flow rate. 
Between stator stagger angle settings of 27.8° and 950° , a 10 per cent 
variation in flow rate requires only an 8° angle change. In the 
regions more distant from the maximum efficiency point, large angle 
changes are required to handle small flow rate changes. 


Figures 31 through 36 also contain plots of 3-stage pressure ratio 


versus referred flow rate for each stator stagger angle. 


S| 


As expected, the measured maximum pressure ratio occurred at a 
referred flow rate between that for maximum efficiency and for surge. 
It was not possible to determine maxima for the calculated data of 
AXCO3 because for each case the largest pressure ratio occurred at 
the flow rate tagged "surge."' The program does not compute perform- 
ance at lower flow rates. The check for surge occurs in subroutine 
CASCAD. Based on cascade data,Howell stated that unacceptable losses 
occur at values of ©/€,, greater than 1.25. This conservative 
criterion is particularly applicable in a design where a wide range 
of flow rates is not required. However for the purpose of evaluating 
pressure ratio for an existing machine, the criterion appears to be 
too restrictive. It would be more advantageous if AXCO3 used the 
criterion to predict surge but then proceeded to complete the calcu- 
lations at that flow rate rather than stopping computation. 

The measured pressure ratio was in general by 0.25 per cent 
greater than predicted by AXCO3, and the measured A'P, ., was about 
11 per cent greater. Total pressure is relatively simple to measure. 
The yaw and pitch errors of the 3-hole probe are very small. Any 
error in total pressure measurement is expected to produce a total 
pressure lower than expected. Since the predictions are lower at 
every flow rate, it is possible that the difference in blading 
between the installed blade profiles and the NACA airfoils, which 
are represented by Liebleins's data, may be the cause of the differ- 
ence as was suggested in the discussion of 3-stage efficiency. A 
summary plot of maximum 3-stage pressure ratio on an expanded scale 
is presented in Fig. 38. The peak value of maximum pressure ratio 


occurred at a stator stagger angle of Suleesb 
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While writing the thesis, an error was discovered in Eq. 1, which 


converts the reading of the mercury barometer Paar in in. of mercury 


to PA in psfa and corrects for the temperature variation of the 


specific gravity of the liquid. The incorrect equation 


P 


= Ve. (0.4928 + Tp, %-9/i0,000)144 (psfa) 


was used. At 70.0°F the constant by which P was multiplied was 


Bar 
71.457 psfa/in. of Hg. Reference 10 lists the specific gravity of 


mercury sy for the temperature range O°F to 150°F as 


sy = 13.638 - 1, 354X10 (T. ) (g/cm?) 


The conversion factor is 


Cy, = Seve (psfa/in. of Hg) 
13,59 


At 70°F the correct conversion factor is 69.55, or the atmospheric 
pressure on which all calculations were based was 3 per cent higher 
than the actual atmospheric pressure. By Eqs. 3, 4 and 7 the cal- 
culated static pressure Fs and air density C, at the permanent 
Pitot-static tube,and the total pressure after the inlet guide vanes 
a ea 3 per cent higher than actual also. By Eqs. 6, 8 and 10 
the calculated velocity Vo? volume flow rate Vee » and referred 
flow rate ah were 0.983 of the actual values. The calculated pres- 


sure ratio P. was 0.998 of the actual value by Eq. 11. By Eqs. 9 


3-s 
and 14 the values of weight flow rate w and 3-stage efficiency 
were calculated to be 1.2 per cent higher than the actual values. 

The test results are presented in non-dimensional form in 
Figs. 39 through 44. Plotted against the flow coefficient ,.,» the 
efficiencies follow the same pattern as when plotted against referred 


flow rate,since Eqs. 10 and 13 differ only by a constant. Figure 45 


is a summary plot of the maximum 3-stage efficiciency for each stator 


5) 


stagger angle versus flow coefficient and is similar to Fig. 37 for 
the same reason. Figure 46 is a summary plot of maximum pressure rise 
coefficient for each stator stagger angle. The pressure rise reaches 
a.maximum at a blade angle of 33° and a flow coefficient of 0.475. The 
curve of corresponding stator stagger angles has an almost constant slope. 
Theoretically, the deflection coefficient Bey represents the 
change in peripheral flow components in a blade row divided by the 
peripheral speed. Figure 47 is a stage velocity triangle showing 
design and off-design conditions in a compressor stage. At the design 


condition 


is) U 


An assumption in this analysis is that the air leaving angle from a 


\N VV 
ls " AVVud _ & Vaud 


blade row does not vary with incidence. Hence the angles &, and (>, 
remain constant. The relation between the off-design deflection 

~! / 
coefficient CL and the off-design flow coefficient ? is fixed by 


similar triangles 


\\ 


and 


or C’ is expected to vary linearly with p- : 

At flow coefficients greater than that for maximum 3-stage 
efficiency,calculated deflection coefficients (Ene varied linearly 
with flow coefficient in Figs. 39 through 44. For most runs the 
graph of ay broke sharply from the linear segment at the maximum 


efficiency point. For each case the curve continued at a reduced 
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slope in the region between the break and prior to surge. Experience 
seems to indicate this identifies a region of incipient stall. As 
the blades become fully stalled, the nature of the flow changes 
drastically. It has been suggested that this change is reflected in 
a sharp increase in the slope of the deflection curve. The author 
was reluctant to operate the compressor in the region of stall for 
extended periods of time because of related vibration and noise prob- 
lems. However a slight upturn of the deflection coefficient curve 
waS apparent at stator stagger angles of 23.8 and 27.8 (Figs. 39 
and 40). 

The deflection coefficient { is a measure of work input. It 
will be necessary to make an approximation for the theoretical power 


Poth’ Equation 11 can be rewritten 
ieee 
ae = | + (ps fa) 
tiga 


Making a first order expansion, Eq. 13 can be rewritten 





’ “7 - ae 
rage =wR, = Jee [ ar a coe | (ft-lbf)/sec 
¥-1 ky; 
1Qa 
w Re T, 
<a eee EEC APS. , = NE oP... (ft-lbf)/sec 
tiga 
Substituting Pith in Eq. 12,the efficiency becomes 
_ Nee aioe 
aS = =e 
aes w 


Rewriting Eq. 15 


T 3 
Vee a EN Pay Uy (f£t~ /sec) 


61 


and introducing Eq. 16 into Eq. 17, there is 


= Va 


iv 
Se,UL" AR Fy 


or the deflection coefficient is a direct function of the torque 
T or work input. 
rq 
Figure 45 is a graph of the ,, which was calculated at max- 
imum 3-stage efficiency for each stator stagger angle. The minimum 
of these deflection coefficients occurs at the maximum efficiency. 


At maximum efficiency the minimum work input was required. 
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SECTION 7 


CONCLUSIONS AND RECOMMENDATIONS 


It has been shown that maximum efficiency can be achieved for 
each flow rate. The computer program AXCO3 can be used to predict 
efficiency within 2 per cent,with an error of 5 per cent in flow 
rate. Pressure ratio may be predicted within 0.25 per cent. 

It is recommended that further measurements be made at increased 
rotor stagger angles. Inlet guide vane angles should be varied to 
simulate the effects of preceding blade rows. Blade shapes with 
higher blade loadings should be tested. 

Scatter in the data may be caused by variations in inlet flow 
and accentuated by oscillations of the pressure gauge. A proposed 


improved inlet duct is described in Appendix E. 


63 





INSTALLATION 


IMPRESSOR 


4YOSS3YdWOD OL S3INOSSI9OV 





2 34xMils 


65 








eae + +oele+t ++44 
tte + a tH 





COMPRESSOR CASING 


FIGURE 3 


S3SITEWSSY LAYNIS ONILNOddNS 
N! GSLNNOW LSAVHS ONY WOHO ‘SS3QV18 YOLON b 3HNOl4 








t ASSEMBLY F ROTOR IN LOWER CASING 





if 


Sa eiNecoceuns SNUNN CET WAYS 


—— 





an 





| 
L 


RE 











g 
Y 





Ai: 


M3IA IWNOILOSS YOSSSudWOD 





Daa 7 















NY 


Ae 
A) of 
iH ol 









— 


ay 
a a a 


Darra 
a 


yd 
N 
a le 














TL VSN? 


TK KA a) | ie 


ZAM 


AEE SBN 3 N N et 
UN Ni iM fa 


9 3yNdlS 


| PR RAPZZZ- 





My a 
IKKE 








69 


S310H AZSAYHNS IWIGWY 4O NOILVOO) 4 ayNdIS 


S3I10H AZAUNS Wilve 3 


SINV dd ABAYUNS JO SNO!LV9D01 ' 
| i Le HI “t. Hy 


Ciimun~ 


oer s¢ (Cie se 
CINVId AZANME WOd NOILDBENO MBIA ———— oe hea el 


IBV ILE ED OL: [| ZZ 


szvez |cosez] sietz | cie'oz aoe 


zat vas out ou 
ws ¢ at % 





cat Lat -——-| 


0€ ,Z1,02 
06 £2.46 


06,25 pri 








70 


YOLOY X3LYOA- 3343 JO 1VL30 8 3YNdI4 





ore | 
mw 27905 
“JOVM aL Jdvm Mass ay 
NOFI ITV WY 10 SILNNIW ES BOIS ggg Sed woLoy vo S/Xv 
AIM ME NOL IIS Jav7yy wosrs FOIs worse ors 
roz7s or3J 


PMO WOILIIC 
NIVI 10 DD 
——~ <340N 








‘VAT OOF” wei ot 


‘wre Ne. awe 


CIVIITY FF Al VION SINL 4O SIKY 
| 4ILON £°o7S 
SOF * aS 727 
' MIIS we AFT! PT 





we Wis ers: 
f/f° 
oee/ 
O04 °2 » 
O99 S 
»® 
oor > 
(4 FL40N FFT} 82% 
Q 
‘IWI0% POAT $ 
DVUNVEL CNV DANMVITN LW MM avy 
ennoarv ; 
gor, IOWTT OL HOw 
= ¥ 00 vore-._¥'" nt pee “tis L00°F AINLIM ALINJOHNI SE 
i aad “SNOILW. 4 
ss! Laer IIS SIINVSI ION heed 
es we OL NULYLNIG 
eee ee ee GIFS INY IN408d JOVIE wOyx'9 
| FVINGIL TOYS Aod IIIS) dtl 
an. GVOAIT NONFT VALI WW SIT MOTT” xz 
vse COVOIL WON SID IGS 


S2VLNMOD FIC WOILIICS® FOS FWSNW) 

IT VIONS YAVIT MW SIWILVIIS FO CLIC 9 
“TIAN 770 he OL AIT 

GNV YCVIT HOLS NSVIS C10 JAowsa Ss 


-5P* OPO -7ITWHVAD 9SOO'* 5492, 
: WINING SILVMEIOO? 
: OL SUNITIOIIV SHOLLIIS Jovrw y¥ 
Sf FLOM’ IIS "FSO AL SBINTOD TIW AINE T 


FOO F SOL 2 "OP-EWN - FIWIS SSINWNONOSY at 
“eFR * O£0°-SBIAWVND WIV 4 
GFIIASIZ IAS FOIMBMIHLOG SSITAMN *S JL0/V 





71 


YOlVWIS X3LNOA-339u4 40 11V1350 ® 3yYNdls 


Qom@ 


mmr ‘o7s xo w23s : 
F995 93S verry wows sav nw 90 SILI EF 
PeNLIM I NOLIIS SOVTE ORs 
















INS Tv F QL FTO SIAL #O STAY 
; ‘340 
w @r Dd w 
Algwyssy av9.sv 
- ind > ae 
t s@-9 3 
AVVO AME :340% 
i 
Sore OO 
een 
oor T 
oor es 
(4 10" 33S} 084 we bbe: 
27 YON I Wie Ware 
POOF S952 
"3WI0Nd 10 ITF 
oot —~\ QI Mah OY DINONIT LY Pw ater awain 735) 
pe SE ; Ganor IW LIINS FOOSE feos 8a! ‘FOOT OL FOVIE wORs 
140°¢ PYALOL WAODAM F/ 


“SN O/LW2 
et(HHIIIS 4£2IV7LNOD JIS 
SIMWAITAL VI/AVLNINO 
MOILIIS TNV WUAGIE POW HES HW 


MOTS Al 





2 oe 
“SN/AA la 
yvasye POL d2 LIS 8) MOE NIAID 
37 OL SQAIIML IT QA LIS ASIN “OH 
leon é $000" NINL IM GY 9 
VA J. , MIia 2 Gav Wlivevd IF RYOHO sian & 
SWOSROONT LOAF JLYNDILIV 7) 9 0864) IIS) Es 010 ‘yOLvIaaVs rags yay 
aaa Ns SS PRET LL Qs® xe yovry 10 SAOSIN OF ULLND 
(F M34 IW » "DUN V0 22d (mm 275) 
0 9; 7 O14 EMAID HNLONFIT VALI ‘NIM TIT MOTI < 


VSWOMLVDITIDPES LIVHLMOD IIS) “MOILIILIA YOS BSW) 
JE PINS WIVIT A! SINILVAEIS OC Sild F 
“FIAO VI" OL SIE 
(NV FOV WOTs NEVI CldW JAOWIE 
"CINSIMPAS SILVMETOOI 
QL DWIAOIDI” SAOLIIS JOVA Y 
“2 510" O4 SUIWVIOD Viv AVIGIT CE 
t 
4% 


“4 


OS S¥M - ¥VMVIS SEIVND0R L 





$00: 519-—h : 
Feo"? of? § O60" -SaF#suUND YI¥ 
GIGI DIGS YOM SEJIAN -$34108 


72 


4qavig 


4YOLVLS X3ILYOA—33y4 


O! 3YNdl4 





os 


IN3SMWLSNray STONY Y3S99VIS YOLOY il 3YuNSls 


ee ars SS -- ——— 
soos 1 
090: 
ESS: 
els" 2 a. ee oe. 
Wl LiNeOLS SIBVTF TIN oF 
INO? C1 - IN 2-9, UI 
¥. Py Ga 
bs aS 
G ) ~ OFSIINIA IP ACW Nia 
HT — - WPI 


YO NOILIIAIO - FLO 
IOY TS /-# i 
a Vig? WRIT C- ; 
YON 0-0 


(IID) GN FOLOI 





a # Py Se 
JOWT@ OL FOVIP SF zt LEO FPO! Lek 


AIVLIS: NO G/F 1-* IQ GFW 
QLITIS” TINS 1-9 


POLO SO 
SINOY 





74 

















FIGURE 12 EXTERNAL STAGGER ANGLE ADJUSTMENT 


= sre 
tt tne , 





FIGURE 13. INLET BELLMOUTH 
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FIGURE 14 
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FIGURE 16 


THREE-HOLE PROBE DETAIL 
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FIG. 1/7 PERMANENT PITOT—STATIC TUBE 






| mm -Ls :. 2 ae | 
FIGURE 18 PERMANENT PITOT ~ TUBE AND 1 HRee— 
HOLE PROBE AHEAD OF FIRST ROTOR 









PRESSURE READOUT INSTRUMENTATION 
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FIGURE 24 TORQUE METER DURING CALIBRATION 
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FIGURE 25 
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FIGURE 26 TRAVERSE OF INLET ODUCT RUN D-7 
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FIGURE 28 TEMPERATURE-ENTROPY DIAGRAM OF 
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FIGURE 30 
3-STAGE EFFICIENCY VERSUS REFERRED 
FLOW RATE PREDICTED BY AXCO 3 
ROTOR STAGGER ANGLE = 43.8° 
VARIOUS STATOR STAGGER ANGLES 
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FIGURE 32 3-STAGE PRESSURE 
RATIO & 3-STAGE EFFICIENCY 
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FIGURE 34 3-STAGE PRESSURE 
RATIO & 3-STAGE EFFICIENCY 
VERSUS REFERRED FLOW RATE 
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FIGURE 358 3-STAGE PRESSURE 
RATIO & 3—STAGE EFFICIENCY 
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APPENDIX A 


CALIBRATION OF INSTRUMENTS 


The measuring devices were calibrated to determine their accuracy 
and precision. 

The Berkeley digital counter is located near a Hewlett-Packard 
digital counter in the control room. At a given flow rate, there 
was no difference between the two instruments. The Hewlett-Packard 
device was recently calibrated with a cesium wave length frequency 
standard in the Standards Laboratory of the Electrical Engineering 
Department of the Naval Postgraduate School. The speed measured is 
considered accurate to 0.5 rpm. 

The Baldwin-Lima-Hamilton torque meter was calibrated statically 
as illustrated in Fig. 24. A bar was attached to the drive end of 
the shaft and was clamped to the stand. A symmetrical lever bar was 
attached to the rotor end of the shaft. The lever arm is 20.0 in. 
Two weight pans were adjusted to 1.44 lb. each with lead shot in a 
plastic bag. The weights themselves were checked on the Toledo 
Precision Scales of the Cascade Laboratory. On the first run, dif- 
ferent readings were obtained while loading and unloading the pans. 
This problem was overcome in subSequent runs by tapping the stand with 
a mallet prior to recording the readdny This simulated the vibration 
of dynamic operation which eliminated any frictional effects. Good 
agreement was obtained during tests where the loading was increased 
and decreased. The possibility of temperature effects on the strain 
gauges was explored in a third test. The torque meter cannot be 


calibrated while the motor is running, but it was calibrated 


Tee 


immediately after shutdown. The shaft was only slightly warmer to 
the touch than ambient temperature. 

The data of the three calibration runs are shown in Table A-l. 
The results are plotted in Fig. A-l. Run A-3 revealed no temperature 
effect. A fourth run was made later for checking purposes and did 
not reveal any discrepancies. The calibration tests established a 
constant for the torque meter of 4.1565 foot-pounds of torque per 
pound reading of the meter. The maximum relative error among the 
runs is 0.6 per cent. The meter may be read to within 0.03 pounds, 
or l per cent of the meter reading. 

The newly acquired low range bourdon tube was calibrated also. 
Run B-l was made by using the mercury manometer board in the control 
room of the Compressor Laboratory. The transsonic turbine test rig 
was used to provide the pressure difference for the higher pressure 
range. The present compressor tests required pressure measurements 
down to 0.6 in. of water. Runs B-3 and B-5 were made at low pres- 
sure difference with a Merriam water micro-manometer. Pressure was 
supplied from a static source. 

The data of the three calibration runs are shown in Table A-2. 
The results are plotted in Fig. A-2. The calculation of velocities 
was set up using pressure difference in wiles Measured pressures 
are therefore very small numbers. An inverse calibration constant 
was used. The bourdon tube has a constant of 240.423 counts per 


2 
lb/ft . The maximum relative error among the calibration runs was 


QO.12 per cent. The meter may be read to within 5 counts. The 


meter's servo system tends to overcorrect. In the manual null mode 


the meter needle tends to wander a similar amount. The five count 
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reading error amounts to 1 per cent of the velocity head reading in 
the inlet duct during flow rate calibration. The error is 0.5 per 
cent of the velocity head reading at the permanent Pitot tube. The 
error is only 0.05 per cent of the measured pressure difference 
across the three stages. 

The Pitot-static tube used to traverse the inlet duct during 
flow rate calibration is a modified Prandtl-type tube. Another 
Pitot-static tube was attached to the traverse shaft and the tubes 
compared. Sufficient difference existed between the two to justify 
further tests. The low speed calibration tunnel in the Cascade 
Laboratory was modified to receive the Pitot-static tube from the 
inlet duct. 

The data of the calibration runs are shown in Table A-3. The 
results are plotted in Fig. A-3. Each point represents the average 
of four readings. Run C-3 indicated that the velocity head of the 
Pitot-static tube from the inlet duct was at most 1.1 per cent less 
than that of the calibration tunnel. Run C-4 indicated that the 
static pressure below atmospheric was at most 1.8 per cent less than 
that of the calibration tunnel. These maximum errors occurred at the 
lower pressure differences. It is probable that the effect recorded 


is in part the reading error due to bourdon tube oscillation. 
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FIGURE A-|! TORQUE METER CALIBRATION 
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: TABLE A-2 ; 
| CALIBRATION OF 0-12 in. Hg BOURDON TUBE | 
= RUN B-1 ba (RUN B-3, | ae RUN B- ae ; 
: a. - Seca | 
PRESSURE METER | PRESSURE i METER : PRESSURE | METER ; 
INCHES OF READING [INCHES OF | READING ‘INCHES OF | READING 
MERCURY COUNTS | WATER COUNTS, WATER COUNTS | 
10.260 178 300 | 13.978 17.409 | 14.519 18,218 | 
8.290 | 139,450 12.459 (15,518 | 12.914 16,119 
7.445 | 125,750| 11.062 13,7840 11.397 14,222 ; 
6.065 102,900 9.853 12,276 ‘=o 12 , 368 : 
5.235 88,200 6.993 , 8,737: 8.471 10, 560 | 
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FIGURE A-2 BOURDON TUBE CALIBRATION 
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TABLE A-3 


CALIBRATION OF INLET PITOT-STATIC TUBE 
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Probe Static Pressure Below 


FIGURE AS INLET TRAVERSE PITOT-STATIC CALIBRATION 





3000 
E 

Run D-3 Probe Crror 
en q Probe — q Tunnel 
5 qa Tunnel 
° 
Oo O 
EB 
z 
= 
8 -Ol 
® 
> 
® 
x=) 
2 
a 

O 1000 2000 3000 
Calibration Tunnel Velocity Head,Counts 
Probe Error 


Pp _ 
SProbe Ps Tunnel 





© Py 
= Tunnel 
3 O 
° 
© 
© 
© 
P= 
a. 
s - Ol 
is 
° 
© 
O 1000 2000 3000 


Calibration Tunnel Static Pressure 
Below Atmosphere, Counts 


121 


APPENDIX B 


FLOW RATE CALIBRATION DETAILS 


Particulars of the flow rate calibration scheme which may be of 
special interest to a restricted number of readers are included in 
this Appendix. Section B-l gives details of the construction and 
testing of the variable length capillary damping devices. Section B-2 
gives a listing of the program ONRFLO, and Section B-3 presents sample 


output of program ONRFLO. 


B-l Variable Length Capillary Damping Devices 


The problem of oscillating pressures at the traversing Pitot- 
static tube in the inlet has been described. Pressure oscillations 
in the tubing were damped with alternating sections of smaller and 
larger diameter flow passages. Figure B-l is a picture of such a 
device. The device consists of three lengths of needle stainless 
steel tubing, of 0.020 in. inside diameter, which are 4 in., 6 in., 
and 8 in. long. They are connected as shown to polyethylene tubing 
by brass "T's.'' By clamping the tubing at various positions it is 
possible to obtain total capillary lengths of 4 in., 6 in., 8 in., or 
18 in. Figure B-2 is a schematic showing the operation of the device. 

Two variable length damping devices were used, one each in the 
total and static pressure lines from the traversing Pitot-static tube 
in the inlet. Two effects of these devices were observed on velocity 
head, the difference between total and static pressure. As the length 


of capillary tube was increased, the fluctuation of the pressure gauge 


was decreased. But as the length of capillary tubing was increased, 


2 


the time to come to a steady reading was also increased. Table B-l 
presents the results of tests of various lengths of capillary tubing. 

Immediately before taking the velocity head reading the traverse 
probe was placed at a new radial position. This task took about 20 
sec. to perform. The best precision for a time to damp of 20 sec 
or less was used. The combination of a total pressure capillary 
length of 18 in. and a static pressure length of 8 in. satisfied this 
requirement. 

Later in the investigation, when the simplified pressure selec- 
tion system was used, it became possible to use both 18 in. lengths 
of tubing and have a time to damp to a steady reading of less than 
20 sec. [It appears that the 12-channel pressure scanners increased 
the time required for the system to stabilize. This was probably due 
to admittance of air from other channels during sequencing of the 
scanner. However as the survey progressed, flow rate measurements 
were made at lower flow rates. As the machine was operated nearer to 
the stalled condition, greater pressure oscillations were observed. 
The meter precision was no better than + 15 counts, and sometimes as 


much as + 50 counts. 


B-2 Program ONRFLO 


The flow rate calculations described in Section 3 were accom- 
plished by program ONRFLO which is included as Table B-2. Since the 
calculations are straightforward and have been described earlier, 


they will not be repeated here. 


The input format is as follows: 


Cageda il: NCASE run number (arbitrary) 
(Iie E7. 2) PBAR reading of barometer (in.) 

TBAR temperature at barometer (°F) 

TO temperature at inlet screen (°F) 
Cards 2, 3, and 4; R the 23 radii at which data were 
(8FLO. 2) taken (17.75, 17.50, ...4.00, 0.00, 


4.00...07250017. 75) @Gine> 
Cards 5 through 50: QI velocity head in inlet duct (counts) 
(4F7.0) PSI static pressure with reference to 
atmospheric in inlet duct (counts) 
QO velocity head at permanent Pitot- 
static tube (counts) 
PSO static pressure with reference to 
atmospheric at permanent Pitot- 


static tube (counts) 


The radii of cards 2, 3 and 4 are repeated internally to account for 
the cross traverse. Cards 5 through 27 introduce measurements made 
at the data points for one traverse, and data for the cross traverse 


are introduced on cards 28 through 50. 


B-3 Output of ONRFLO 


Table B-3 is a sample of the output of ONRFLO. It is for the 
final calibration, run 14. The first data block is a list of input 
data. For each point on the two traverses, QI, PSI, QO, and PSO are 
listed. The second data block is a list of computed velocities VI, 


VO, and VIAD. The third block shows how the "no boundary layer" 


124 


flow rate, VFI, is computed by summing the product of VA and A. 
The fourth block shows how the flow rate deficiency, VFD, is computed. 
It lists the inputs to and outputs from SUBROUTINE QTFE, which per- 
forms the integration by the trapezoidal rule. The final statements 
list the final computed flow rate VFC and the final calibration 


constant CCC. 
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FIGURE B-i VARIABLE LENGTH CAPILLARY DAMPING DEVICE 
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Listing of Program ONRFLO (cont. ) 
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Listing of Program ONRFLO (cont. ) 
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Sample Output of Program ONRFLO-Input for Run 14 
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Sample Output of Program ONRFLO-Velocities for Run 14 
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APPENDIX C 


PERFORMANCE MEASUREMENTS 


It was decided to take flow measurements at the mean streamline, 
using the radial survey holes. The assumption was made that measure- 
ments made at this point accurately represented the flow in the 
entire flow annulus. The truth of the assumption is verified in 
Section C-l. Section C-2 presents details concerning the data re- 


duction program ONRETA. 


C-1 Verification of Assumptions 


The assumption that measurements made at the mean streamline 
using the radial survey holes accurately represent the flow in the 
entire flow annulus was verified in three ways. The total pressure 
measured at the mean radius was shown to be nearly equal to the 
integrated average of total pressure at each radius. It was demon- 
strated that the total pressure measured at the peripheral location 
of a radial survey hole was nearly equal to the integrated average of 
total pressure at each peripheral location in a blade channel. The 
flow was shown to be nearly axisymmetric. 

Concerning the radial variation of total pressure, the existence 
of a boundary layer at the outer radius in the inlet has already been 
shown. It was further shown that a boundary layer existed ahead of 
the first stage and behind the last stage. With the stator stagger 
angle at 27.8 and the flow rate at about 400 Fae ec. a radial 
traverse was made at the radial survey hole labeled 37°27' 30" in 


survey plane 8 of Fig. 7. Data points were established each 0.3 in. 
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from a radius of 11.1 in. to 18.0 in. After recording the barometric 


pressure P and ambient temperature se the following measurements 


Bar 
were taken: 
Pid - total pressure with reference to atmospheric after third 
Stator (counts) 
Pod - static pressure with reference to atmospheric after third 


stator (counts) 


qq 7 velocity head after the third stator (counts) 


The total pressure at the inner radius could not be measured but was 
assumed to be the static pressure at 11.1 in. With the relations 


listed on p. 35 these measuring data establish the following 


quantities; 
P= Pao 071.467) (ps fa) (1) 
T= Ty + 4547 CR) (2) 
Ped 
ead Valea 240.423 (psfa) Sy 
Psd 
reaeie 379553 a si 


With the assumption that the total temperature after the third stator 


ways is nearly equal to the inlet total temperature Teo 


p ’ 
e, ee ae ST eda 518.7 sta (5.47 x10") an sec’/ft”) (4) 
(14,696)(144) Ty Teo 
, Id 
bdi 240.423 ee) 


136 


= dh (ft/sec) (6) 


Ra 


The average total pressure was calculated from the integral expression 


{8.0 
ae: d 
a i 7 Eig % 
> 10,8 
re = 
18,0 


yell V, a ie 
10.8 


Both integrations were performed by the trapezoidal rule. Figure C-l 
is a graphical representation of the integration. The integrated 
average total pressure ae was 2184.20 Lb/£t? while the total pressure 
measured at the mean streamline was 2184.34 yet. The relative 
error between the two pressures is 0.006 per cent. A second radial 
survey was made at the radial survey hole labeled 149°57'30" in the 
plane called S. P. 8 on Fig. 7. The integrated average pressure for 
this survey was 2183.92 lee The total pressure measured at the 
mean radial position was 2184. 32 infers The relative error between 
the two pressures was 0.018 per cent. Since the pressure measured 

at the mean radius was nearly equal to the integrated average of pres- 
sure at each radial location, mean streamline values were used to mea- 
sure compressor performance. 

The traverse carriage was used for peripheral surveys. Since it 
would be necessary to remove the first stage rotor to survey in the axial 
plane ahead of the first stage, the investigation was confined to the 
plane behind the third stage stator where measurements could be made simply. 
The locations of the radial survey holes with reference to the indicated 


meridional angle of the carriage are shown in Fig. 23. 
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There are 32 stator blades. Each blade channel is 11.25° in peripheral 
direction. With the stator stagger angle at ieGmae peripheral survey 
was made at the mean radius, 14.4 in. The average reading over the 
entire 15° of carriage travel was 9256 counts. The reading taken at 
the radial survey hole marked 37°27'30" on Fig. 7 was 9246 counts, a 
relative error of 0.11 per cent. The average over the first 11.25° of 
carriage travel, representing one blade channel, was 9218 counts, a rel- 
ative error of 0.30 per cent. Figure C-2 illustrates the results of 
this mean radius survey. A second survey was conducted near the hub 
at 2 eydtes of 11.2 in. The average pressure reading of the traverse 
was 9532.47 counts. The pressure measurement at the radial survey hole 
was 9569.23 counts. A relative error of 0.36 per cent existed between 
the two pressures. On the basis of the surveys made behind the third 
Stage stator the measurements taken at the radial survey taps were 
accepted as accurately representing the flow in the entire blade 
channel. 

An investigation of the axisymmetry of the flow was conducted. 
The results are tabulated in Table C-l. The total pressure behind the 
inlet guide Ce, Ue eae measured at eight different radial survey 


holes in plane S. P. 2 of Fig. 7. The average reading of P was 


tig 
578.1 counts. The nearest reading to the average was 573.9 counts, 

taken at the hole labeled 33°45'. This hole was used for all readings 
of Ye The relative error was 0.73 per cent. The average of read- 


ings taken behind the third stage stator las was 10, 821.5 counts. 


d 


The nearest reading to the average was 10, 810 counts, taken at the 


hole marked 37°27'30" on Fig. 7. This hole was used for all readings 


of Pea’ The relative error was 0.10 per cent. Within the measurement 
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capability of the instrument, the flow was considered axisymmetric 


and the most representative readings have been used. 


C-2 Program ONRETA 


The data reduction described in Section 4 was accomplished by 


Program ONRETA (Table C-2). 


Since the calculations are straight- 


forward and have been described in the body of the thesis, they will 


not be repeated here. 


The input format is as follows: 


Card 1: NCASE 
(218 ,4F8.2) NRUN 
PBAR 
TBAR 
RSTAG 
SSTAG 
Cara 2: TO 
(F8.1,4F8.0 QO 
F8.1,F8.2) 


PSO 


PTIG 


DELP 


TACH 
TRAW 


Cards 3 through (NRUN - 1): 


run number (arbitrary) 

number of data points per run 
reading of barometer (in. of Hg) 
temperature of barometer (°F) 
rotor stagger angle (43.8°) 
Stator stagger angle (various a) 
temperature at inlet screen (°F) 
velocity head at permanent Pitot- 
static tube (counts) 

static pressure with reference to 
atmospheric at permanent Pitot- 


static tube (counts) 


total pressure with reference to atmo- 
spheric ahead of first rotor (counts) 


pressure rise across the 3 stages 
(counts) 


reading of the speed counter (rpm) 


reading of the torque meter (1b) 


repeat of Card 2 for remaining data points 
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The first five entries on cards 2 through (NRUN - 1) are the arithmetic 
average of four readings taken two minutes apart. The last two entries 
are the average of two readings taken five minutes apart. 

The first section of program output is a printout of the input 
data. The second section is the calculated data for the run. The 
third section is the calculated nondimensional coefficients. 


Tables C-3 through C-13 present the results of runs 1 through 11. 
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Listing of Program ONRETA 
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Output of Program ONRETA, Run 1, Stator Stagger Angle = 27.8 
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Table C-6, Output of Program ONRETA, Run 4, Stator Stagger Angle = 23.8. 
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Table C-11. Output of Program ONRETA, Run 9, Stator Stagger Angle = 35.8 
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Output of Program ONRETA, Run ll, Stator Stagger Angle = 44,3° 
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APPENDIX D 


DETAILS OF PREDICTION PROGRAM 


This Appendix contains details concerning AXCO3 which may be of 
interest to a limited number of readers. Table D-1l is a listing 
of the main program and the seven subroutines. Section D-l is a 
description of the method of inserting data. Table D-2 is a sample 
of program output. The rewrite of input data, a sample of interstage 
data,and a sample summary of one case are included. Tables D-3 
through D-16 are the summary printouts of the cases which were compared 
to measured performance of the machine. Table D-17 is a printout of 
program CHECK which verified the errors in the polynomials of sub- 


routines THEORY and CASCAD. 


D-1 Input Data 


There is no limit to the number of cases which may be considered, 
since data are read in prior to execution for each case, and output 
printed prior to considering the next case. The executive routine 
reads input data and prints it before computation begins. The input 


data consists of thirteen cards per case. 


Card l: NUM case number (arbitrary) 
(l73 GAM ratio of specific heats 
4F7. 3) RN referred speed (xpm/V°R) 
WREM initial referred flow rate 
(1bm\f°rR/sec psia) 
DW referred flow rate increment 


(1bm>\/ °R/sec psia) 
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Card 2: 
Gi? ; 


SF7. 3) 


Card 3: 
(3F7.3, 
457.1, 


3F7. 3) 


Card 4; 
(10F7. 3) 
Card 5; 
(10F7. 3) 
Card 6; 
(10F7. 3) 
Card 7: 


(10F7. 3) 


LMAX 


ZETAL 


ZIGV 


EGV 
ETD 
EEX 
DIN 
DAV 


DOU 


AIGV 
AGV 
ADI 
ALE 
SMAX 
SMI 


AAA 


BH 


BLOCK 


WDO 


number of stages (3) 

loss coefficient of inlet duct 

(assumed 0.0) 

loss coefficient of inlet guide vanes 
(assumed 0.05) 

efficiency of exit guide vanes (assumed 0.85) 
efficiency of diffuser (assumed 1.0) 
efficiency of exit (assumed 1.0) 

inlet pipe diameter, in. 

average blade diameter, in. 

outlet pipe diameter, in. 

annulus area ahead of IGV, ee 

annulus area after IGV, ae 

annulus area after EGV, inne 

diffuser exit area, see 

inlet guide vane exit angle, degrees 
maximum limit, Ca) ae set by Fig. 29 
minimum limit, (tal a duhenie » set by Fig. 29 
an array of annulus areas after rotor one, 
stator one,....stator LMAX, ine 

an array of blade heights of rotor one.... 
stator LMAX, in. 


an array of blockage factors for rotor one.... 


stator LMAX (assumed 1.0) 


an array of work done factors for stage one.... 


stage LMAX (assumed 1.0) 
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Card 8: 
(Rotor 
data) 


(9F7.3) 


Card 9: 
(Stator 
data) 


QOF7. 3) 


BLDR 


CAMR 


STAGR 


SHR 


THR 


CHR 


DIAR 


DER 


CDR 


BLDS 


CAMS 


STAGS 


SHS 


THS 


CHS 


DIAS 


DES 


CDS 


number of blades in the row 

camber angle, degrees 

Stagger angle, degrees 

shape correction factor, (1.0) 

thickness to chord ratio 

blade chord length, in. 

mean blade diameter, in. 

tip clearance, in 

minimum profile drag coefficient (assumed 
various) 

number of blades in the row 

camber angle, degrees 

stagger angle, degrees 

shape correction factor (1.0) 

thickness to chord ratio 

blade chord length, in. 

blade mean diameter, in. 

tip clearance, in. 

minimum profile drag coefficien_ (assumed 


various) 


Cards 10 and 12 are duplicates of 8 


Cards 11 and 13 are duplicates of 9 
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OVERALL DATA 
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Output of Program AXCO3, Case 32, Cp; 
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Table D-11. Output of Program AXCO3, Case 33, Cys, 0.008 
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Table D-12. Output of Program AXCO3, Case 41, Comin 2:000 
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Output of Program AXCO3, Case 42, Cys, 0.006 
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Output of Program AXCO3, Case 43, C_. 


Table D-14. 
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65eUUNUU00V000 
7c. JOUVU0UU9NGO 


cs UVEGrCES 


we 

2 eVJUNDVUD0UU0 
Lu OQV0VU00V0V0U 
L5.-QOQ0VU00V0000 
¢ }eCOUDUOVVNO 
é2.UGU00000VU0 
$3.2 UQOUCQUCUN 
35.-QVUUVUCO0VU0 
1 eVGVUUNV0VE 
+5eur0CU0U0N0 
»0.0UVIU000U0 
»o9 YVIICUUGUUYU 
o0eVIVVIVUNIUNYU 
639 VJUNUUVV0V0N 
fo-GCUuuQ0U00 


ZLKG CAMBER INCIDENCE VERSUS ANGLE 
ACTUAL INCIDENCE CALCULATEO INCIDENCE 


U.0 

Ue 3299999 833 
0.6799999475 
09899999499 
1. 3099994059 
1-6199998856 
1.9299993515 
202299995422 
225499992371 
2.8299999237 
3-1199996556 
323899993896 
3-6199996356 
325299999237 
4.0099992752 


SLOPE FACTOR N VERSUS ANGLE 


ACTUAL FACTOR 


-0.0389999971 
—-0.0499959970 
-J.0609999970 
~J20739999413 
-0.057999S9995 
—-J. 1039999 7/23 
-UV- 1209999919 
— 321499999725 
-J- 1619999409 
-Geldd999G895 
~J-2119999528 
—U.2409999967 
-J-216999S504 
—~VU.- 3209999500 
-J.37599996 71 


OEVLATIUN VERSUS ANGLE 


eu 

0.0599999987 
0-129999S952 
U.20999S9T 55 
0. 29999599523 
Je 394999959762 
Je-508S9999905 
Ve bo 99999952 
UO. 7539999905 
J- 3999995 762 
1.0599994059 
1.250000UU0 UU 
1.4799995422 
1. 7500000VU00 
2-U69999¢948 


Table D- 17. 


SOLIDITY = 0.8 


-0.0009329878 
0.3080989718 
0.6169836521 
0-9257200956 
1.2343072891 
125427494049 
1.8510446548 
221591939926 
224671907425 
227750444412 
3.0827503204 
33903064728 
3.6977157593 
4.-0049800873 
423120946884 


CALCULATED FACTOR 


—-0 .0343636572 
-0.0483931080 
-0.0621372499 
-0.0760737062 
-0.09068071 84 
—-0.1064355373 
-~0.-1238162518 
—-0.1433003545 
-0-1653655767 
-0.1904897094 
~0.2191501260 
-0.2518250346 
-0.2889920473 
~0 «3311289549 
—0.3787130713 


ACivLaAl UcVIATIC’ CALCULATED DEVIATION 


—-0.0014475845 
0.0814747810 
0.1593248316 
0.2358062267 
0-3145923018 
0.3993784785 
0.4938524365 
0.601 7042994 
0. 7266232967 
0.8722926378 
1-0424118042 
1.2406597137 
1.4707336426 
1.7363204956 
2 -0411033630 
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ACTUAL - 


0-0009329878 
0.0219010115 
0-0630162954% 
0.0642798543 
0-0756921768 
0.0772504807 
0.0789546967 
0.0708055496 
0.0828084946 
020549554825 
0.0372495651 
~0.200030 70831 
—-0.0777158737 
-0.1749801636 
-0.3020954132 


ACTUAL — CALCULATED 


-~0.0046363398 
-0-0016068891 
0-0011372529 
0.0020737648 
0-.0026807189 
0.0024355650 
0.0026162599 
0.0023003817 
0.0033656359 
0-0044897199 
0.0071501 732 
0.0108250380 
0-0119920969 
0.0101289749 
0.0027130842 


ACTUAL — DEVIATION, 


0-0014475845 
—-0.0214747824% 
-0.0293298364 
—-0-0258062482 
-0.0145923495 

020006214976 

0.01614 75539 

0.0282956958 

00333766937 

0.0277073383 

0-0175876617 

0.0093402 863 

0.0092658997 

0.0136795044 

0.0288963318 


Program CHECK,Listing and Output (cont. ) 


CALCUL ATEO 


PERCENT ERROR 


0-0932987332 
66366682053 
922671022415 
6 «4929141998 
5.7780303955 
4-7685480118 
4-0909175873 
321751365662 
322473926544 
1.9418897629 
1-1938962936 
-0.00905 84978 
—221468467712 
~4.5686721802 
=7 65335502625 


PERCENT ERROR 


11.8880500793 

322137775421 
~1-8643474579 
~228023862839 
-—3 0462713242 
—-2-3418893814 
—2.3274869919 
~1.-6314764023 
~2.0775537491 
~264138278961 
—3.3727235794 
~4.4917144775 
=4 «3292791512 
- 321554431915 
-0-7215649486 


PERCENT ERROR 


0.1447584033 
~35.7912902832 
~22.5614013672 
~12.28868961 33 

—4-.8641147614 

0.1553744078 

3- 1661853790 

44913787842 

423916702271 

3.0785923004 

126592140198 

07472229004 

06260744929 

0.7816859484 

1.3959579468 


ANGLE, CEGRECS 


0.0 

2 eVUVD000000 
1¢.0000000000 
15.00V0000000 
2G-0000000000 
2520000000000 
3UeOUDUCVLVUN 
35.-VUV00V00000 
40- 0000000000 
45.0000000000 
50-0000000000 
55 -UV0VVU00000 
60.000D00U000 
65e¢00Q0ULI00U 
70.0000000000 


ANGLts ULLLKELS 


0.0 
2-VQUVU0000U 
L0.0000000U00 
15.0000000000 
20 -GC00VU000V000 
2>.-VJ0V0C0000 
30.0000000000 
33-0V00000VU00 
4uU-00Q000000U00 
45.UU00000000 
5G -QUVD000000 
25-VOU0U0U0I00U 
b606«V00U00000J0 
t52QVQV0000000 
f0.0J900000U0C0 


SLUPE FACTOR M VERSUS ANGLE 


ACTUAL FACTOR ACTUAL - CALCULATED 


0.3099999428 
0. 3119999766 
0. 3149999976 
0. 3169999719 
YQ. 3209999800 
0-3239999413 
0.3279999495 
0. 3329999447 
0.3379999995 
0.3439999819 
0.3499999642 
0. 3569999933 
0. 3649999499 
0.375Q0000000 
0-3879999518 


CALCULATED FACTOR 


0.3104162812 
0.3125761747 
0-3147908449 
0.3171648979 
03198031187 
0. 3228106499 
0.3262923956 
0-3303530812 
0.3350979686 
0-3406320214 
0.3470600843 
0.3544869423 
0-3630179167 
003727575541 
0.3838108182 


-0-0004163384 
-0-0005761981 


0-0002091527 


-0-0001649261 


0.0011968613 
0.0011892915 
0-0017075539 
0-0026468635 
0.0029020309 
0.0033679605 
0.0029398799 
0.0025130510 
0.0019820333 
0-0022424459 
0.0041891336 


SLOPE FACTOR 01/00 VERSUS ANGLE/ 


ACTUAL FACTOR 


0.0769999623 
0.0789999962 
0.0809999704 
0.082999S9447 
0.085999965/7 
0.0899999738 
0.09399¢99819 
0.1019999981 
0-110999941%& 
U.1209999919 
0.1319999095 
V2~1499999 162 
0.1719999909 
0.-1979999542 
We 2269599981 


Table D-17. 


CALCULATED FACTOR 


0.075421 8698 
0.078034 7586 
0.0802936554 
0.0826053023 
0-0853770375 
0-0890154839 
0.0939274430 
0.1005203128 
0.1092003584% 
0.1203768849 
0.1344547272 
0-1518397331 
021729413271 
0.1981655955 
U.-2279177904 


195 


- CALCULATED 


0.0015780926 
0.0009652376 
0-0007063150 
020003946424 
0-0006229281 
0-0009844899 
0.0000725389 
0.0014 796853 
0.0017995834 
0-0006231070 


—-0.0024547577 
-0.0018397570 
-0.0009413362 
-J0.000 1656413 
-0.0009177923 


Program CHECK,Listing and Output (cont.) 


PEKCENT ERKOR 


-0.1343027353 
~0-1846788526 


0.06639 76669 


-0.0520271473 


0-3728539944 
0.3670653105 
0-5205957294 
0.7948539853 
0-8585889935 
029790582657 
0.8399657607 
0.7039358616 
0.5430228710 
0-5979855657 
1.07967281 34 


PERCENT ERROR 


2-049471865527 
122218189240 
0.8719941373 
004754729667 
0.7243353128 
1.0938777924 
0.0771690011 
1.4506711960 
1.6212463379 
0.5149644613 


—-1.-8596649170 
-1.2265043259 
—-0.5472884774 
-0.0836572051 
—0-4043137431 


APPENDIX E 


PROPOSED DESIGN OF AN IMPROVED INLET DUCT 


The present inlet duct is deficient because of the shape of its 
bellmouth, because of the small velocities that occur in it during 
operation which make accurate determination of volume flow rates 
impossible, and because it has no provision to eliminate effects of 
atmospheric gusts or winds. 

The entrance bellmouth should have more gradual changes in wall 
curvature to eliminate local separations at the entrance to the 
cylindrical duct. It was stated earlier that the velocity head in 
this duct with 36 in. diameter is about 0.6 in. of water, a value too 
small for accurate measurements. To avoid large fluctuations of the 
read-out of the Texas Instruments pressure gauge, for small variation 
in flow rate, the velocity head should be increased to about 6 in. 
of water or 33 ey et", at a volume flow rate of 360 ft” ee giving 
a response at the instrument of about 8000 counts. For an assumed 
incompressible flow with a mass density Q 2 of 0.002378 lb mT, 
ae and a velocity head qi of 33 ‘n/a the velocity V, in the 


a8 


measuring plane of the inlet duct would have to be 


A Hit 
Vee | ——- (ft/sec) (6) 


L a 


Hence the required flow area is 360/166 = 2.16 ae corresponding 





to a circular area with 19.9 in. diameter. Because of the small 
pressure rise of the compressor it is necessary to convert the 


velocity head at the measuring plane into static pressure rise by 


arranging an optimum diffuser between this section and the inlet 
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pipe to the compressor that has a diameter of 36 in. Test data of 
Reference 11 show that the included diffuser angle should be 8° for 
the present application. Larger diffuser angles produce separation 
with associated flow instabilities, and smaller angles require ducts 
with increased lengths where the boundary layer growth is excessive 
and reduces pressure recovery. For a diffuser angle of 8° the length 


L of the diffuser is 


| = ae 14.9 = WIG tin, 


2 tan 4 


Figure E-l shows the design of the proposed inlet duct. A well 
designed bellmouth from a diameter of 48 in. to the throat of 19.9 
in. requires an axial length of about 26 in., giving a total length 
of the duct of 116 + 26 = 142 in. Ahead of the bellmouth inlet a 
structure will be attached to support several layers of fine mesh 
screen to eliminate flow disturbances by gusts in the surrounding 
atmosphere. The whole duct is mounted on a trolley which can be 
rolled onto the apron outside the test cell to permit simple attach- 
ment to and removal from the presently installed cylindrical inlet 
duct. Jt is recommended also that an additional honeycomb flow 


Straightener be installed at the entrance of this duct. 


The bellmouth and the diffuser could be molded out of plastics, 
reinforced by fiberglass, by using a wooden template which is split 


at the smallest diameter. 
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